Integrative analysis of DNA replication origins and ORC-/MCM-binding sites in human cells reveals a lack of overlap.

Based on experimentally determined average inter-origin distances of ~100 kb, DNA replication initiates from ~50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the origin recognition complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and five ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ~7.5 million union origins identified by all datasets, only 0.27% (20,250 shared origins) were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques, suggesting extensive variability in origin usage and identification. Also, 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF-binding sites, G-quadruplex sites, and activating histone marks, these overlaps are comparable or less than that of known transcription start sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ~13,000 reproducible ORC-binding sites in human cancer cells, and only 4.5% were within 1 kb of the ~11,000 union MCM2-7-binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, Saccharomyces cerevisiae. Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity.

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

DARDN: A Deep-Learning Approach for CTCF Binding Sequence Classification and Oncogenic Regulatory Feature Discovery.

Characterization of gene regulatory mechanisms in cancer is a key task in cancer genomics. CCCTC-binding factor (CTCF), a DNA binding protein, exhibits specific binding patterns in the genome of cancer cells and has a non-canonical function to facilitate oncogenic transcription programs by cooperating with transcription factors bound at flanking distal regions. Identification of DNA sequence features from a broad genomic region that distinguish cancer-specific CTCF binding sites from regular CTCF binding sites can help find oncogenic transcription factors in a cancer type. However, the presence of long DNA sequences without localization information makes it difficult to perform conventional motif analysis. Here, we present DNAResDualNet (DARDN), a computational method that utilizes convolutional neural networks (CNNs) for predicting cancer-specific CTCF binding sites from long DNA sequences and employs DeepLIFT, a method for interpretability of deep learning models that explains the model's output in terms of the contributions of its input features. The method is used for identifying DNA sequence features associated with cancer-specific CTCF binding. Evaluation on DNA sequences associated with CTCF binding sites in T-cell acute lymphoblastic leukemia (T-ALL) and other cancer types demonstrates DARDN's ability in classifying DNA sequences surrounding cancer-specific CTCF binding from control constitutive CTCF binding and identifying sequence motifs for transcription factors potentially active in each specific cancer type. We identify potential oncogenic transcription factors in T-ALL, acute myeloid leukemia (AML), breast cancer (BRCA), colorectal cancer (CRC), lung adenocarcinoma (LUAD), and prostate cancer (PRAD). Our work demonstrates the power of advanced machine learning and feature discovery approach in finding biologically meaningful information from complex high-throughput sequencing data.

The transcriptional cofactor Tle3 reciprocally controls effector and central memory CD8+ T cell fates.

Antigen-experienced CD8+ T cells form effector and central memory T cells (TEM and TCM cells, respectively); however, the mechanism(s) controlling their lineage plasticity remains incompletely understood. Here we show that the transcription cofactor Tle3 critically regulates TEM and TCM cell fates and lineage stability through dynamic redistribution in antigen-responding CD8+ T cell genome. Genetic ablation of Tle3 promoted CD8+ TCM cell formation at the expense of CD8+ TEM cells. Lineage tracing showed that Tle3-deficient CD8+ TEM cells underwent accelerated conversion into CD8+ TCM cells while retaining robust recall capacity. Tle3 acted as a coactivator for Tbet to increase chromatin opening at CD8+ TEM cell-characteristic sites and to activate CD8+ TEM cell signature gene transcription, while engaging Runx3 and Tcf1 to limit CD8+ TCM cell-characteristic molecular features. Thus, Tle3 integrated functions of multiple transcription factors to guard lineage fidelity of CD8+ TEM cells, and manipulation of Tle3 activity could favor CD8+ TCM cell production.

Proximal immune-epithelial progenitor interactions drive chronic tissue sequelae post COVID-19.

The long-term physiological consequences of SARS-CoV-2, termed Post-Acute Sequelae of COVID-19 (PASC), are rapidly evolving into a major public health concern. The underlying cellular and molecular etiology remain poorly defined but growing evidence links PASC to abnormal immune responses and/or poor organ recovery post-infection. Yet, the precise mechanisms driving non-resolving inflammation and impaired tissue repair in the context of PASC remain unclear. With insights from three independent clinical cohorts of PASC patients with abnormal lung function and/or viral infection-mediated pulmonary fibrosis, we established a clinically relevant mouse model of post-viral lung sequelae to investigate the pathophysiology of respiratory PASC. By employing a combination of spatial transcriptomics and imaging, we identified dysregulated proximal interactions between immune cells and epithelial progenitors unique to the fibroproliferation in respiratory PASC but not acute COVID-19 or idiopathic pulmonary fibrosis (IPF). Specifically, we found a central role for lung-resident CD8+ T cell-macrophage interactions in maintaining Krt8hi transitional and ectopic Krt5+ basal cell progenitors, thus impairing alveolar regeneration and driving fibrotic sequelae after acute viral pneumonia. Mechanistically, CD8+ T cell derived IFN-γ and TNF stimulated lung macrophages to chronically release IL-1ß, resulting in the abnormal accumulation of dysplastic epithelial progenitors and fibrosis. Notably, therapeutic neutralization of IFN-γ and TNF, or IL-1ß after the resolution of acute infection resulted in markedly improved alveolar regeneration and restoration of pulmonary function. Together, our findings implicate a dysregulated immune-epithelial progenitor niche in driving respiratory PASC. Moreover, in contrast to other approaches requiring early intervention, we highlight therapeutic strategies to rescue fibrotic disease in the aftermath of respiratory viral infections, addressing the current unmet need in the clinical management of PASC and post-viral disease.

Integrative single-cell meta-analysis reveals disease-relevant vascular cell states and markers in human atherosclerosis.

Coronary artery disease (CAD) is characterized by atherosclerotic plaque formation in the arterial wall. CAD progression involves complex interactions and phenotypic plasticity among vascular and immune cell lineages. Single-cell RNA-seq (scRNA-seq) studies have highlighted lineage-specific transcriptomic signatures, but human cell phenotypes remain controversial. Here, we perform an integrated meta-analysis of 22 scRNA-seq libraries to generate a comprehensive map of human atherosclerosis with 118,578 cells. Besides characterizing granular cell-type diversity and communication, we leverage this atlas to provide insights into smooth muscle cell (SMC) modulation. We integrate genome-wide association study data and uncover a critical role for modulated SMC phenotypes in CAD, myocardial infarction, and coronary calcification. Finally, we identify fibromyocyte/fibrochondrogenic SMC markers (LTBP1 and CRTAC1) as proxies of atherosclerosis progression and validate these through omics and spatial imaging analyses. Altogether, we create a unified atlas of human atherosclerosis informing cell state-specific mechanistic and translational studies of cardiovascular diseases.

Proximal immune-epithelial progenitor interactions drive chronic tissue sequelae post COVID-19.

The long-term physiological consequences of SARS-CoV-2, termed Post-Acute Sequelae of COVID-19 (PASC), are rapidly evolving into a major public health concern. The underlying cellular and molecular etiology remain poorly defined but growing evidence links PASC to abnormal immune responses and/or poor organ recovery post-infection. Yet, the precise mechanisms driving non-resolving inflammation and impaired tissue repair in the context of PASC remain unclear. With insights from three independent clinical cohorts of PASC patients with abnormal lung function and/or viral infection-mediated pulmonary fibrosis, we established a clinically relevant mouse model of post-viral lung sequelae to investigate the pathophysiology of respiratory PASC. By employing a combination of spatial transcriptomics and imaging, we identified dysregulated proximal interactions between immune cells and epithelial progenitors unique to the fibroproliferation in respiratory PASC but not acute COVID-19 or idiopathic pulmonary fibrosis (IPF). Specifically, we found a central role for lung-resident CD8+ T cell-macrophage interactions in maintaining Krt8hi transitional and ectopic Krt5+ basal cell progenitors, thus impairing alveolar regeneration and driving fibrotic sequelae after acute viral pneumonia. Mechanistically, CD8+ T cell derived IFN-γ and TNF stimulated lung macrophages to chronically release IL-1ß, resulting in the abnormal accumulation of dysplastic epithelial progenitors and fibrosis. Notably, therapeutic neutralization of IFN-γ and TNF, or IL-1ß after the resolution of acute infection resulted in markedly improved alveolar regeneration and restoration of pulmonary function. Together, our findings implicate a dysregulated immune-epithelial progenitor niche in driving respiratory PASC. Moreover, in contrast to other approaches requiring early intervention, we highlight therapeutic strategies to rescue fibrotic disease in the aftermath of respiratory viral infections, addressing the current unmet need in the clinical management of PASC and post-viral disease.

Integrative analysis of DNA replication origins and ORC/MCM binding sites in human cells reveals a lack of overlap.

Based on experimentally determined average inter-origin distances of ∼100 kb, DNA replication initiates from ∼50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the Origin Recognition Complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and 5 ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ∼7.5 million union origins identified by all datasets, only 0.27% were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques (20,250 shared origins), suggesting extensive variability in origin usage and identification. 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF binding sites, G-quadruplex sites and activating histone marks, these overlaps are comparable or less than that of known Transcription Start Sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ∼13,000 reproducible ORC binding sites in human cancer cells, and only 4.5% were within 1 kb of the ∼11,000 union MCM2-7 binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, S. cerevisiae . Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

Single cell RNA sequencing unravels mechanisms underlying senescence-like phenotypes of alveolar macrophages.

Alveolar macrophages (AMs) are resident innate immune cells that play vital roles in maintaining lung physiological functions. However, the effects of aging on their dynamics, heterogeneity, and transcriptional profiles remain to be fully elucidated. Through single cell RNA sequencing (scRNA-seq), we identified CBFß as an indispensable transcription factor that ensures AM self-renewal. Intriguingly, despite transcriptome similarities of proliferating cells, AMs from aged mice exhibited reduced embryonic stem cell-like features. Aged AMs also displayed compromised DNA repair abilities, potentially leading to obstructed cell cycle progression and an elevation of senescence markers. Consistently, AMs from aged mice exhibited impaired self-renewal ability and reduced sensitivity to GM-CSF. Decreased CBFß was observed in the cytosol of AMs from aged mice. Similar senescence-like phenotypes were also found in human AMs. Taken together, these findings suggest that AMs in aged hosts demonstrate senescence-like phenotypes, potentially facilitated by the abrogated CBF ß activity.

Fob1-dependent condensin recruitment and loop extrusion on yeast chromosome III.

Despite recent advances in single-molecule and structural analysis of condensin activity in vitro, mechanisms of functional condensin loading and loop extrusion that lead to specific chromosomal organization remain unclear. In Saccharomyces cerevisiae, the most prominent condensin loading site is the rDNA locus on chromosome XII, but its repetitiveness deters rigorous analysis of individual genes. An equally prominent non-rDNA condensin site is located on chromosome III (chrIII). It lies in the promoter of a putative non-coding RNA gene called RDT1, which is in a segment of the recombination enhancer (RE) that dictates MATa-specific chrIII organization. Here, we unexpectedly find that condensin is recruited to the RDT1 promoter in MATa cells through hierarchical interactions with Fob1, Tof2, and cohibin (Lrs4/Csm1), a set of nucleolar factors that also recruit condensin to the rDNA. Fob1 directly binds to this locus in vitro, while its binding in vivo depends on an adjacent Mcm1/α2 binding site that provides MATa cell specificity. We also uncover evidence for condensin-driven loop extrusion anchored by Fob1 and cohibin at RDT1 that unidirectionally extends toward MATa on the right arm of chrIII, supporting donor preference during mating-type switching. S. cerevisiae chrIII therefore provides a new platform for the study of programmed condensin-mediated chromosome conformation.

Interpretable meta-learning of multi-omics data for survival analysis and pathway enrichment.

MOTIVATION: Despite the success of recent machine learning algorithms' applications to survival analysis, their black-box nature hinders interpretability, which is arguably the most important aspect. Similarly, multi-omics data integration for survival analysis is often constrained by the underlying relationships and correlations that are rarely well understood. The goal of this work is to alleviate the interpretability problem in machine learning approaches for survival analysis and also demonstrate how multi-omics data integration improves survival analysis and pathway enrichment. We use meta-learning, a machine-learning algorithm that is trained on a variety of related datasets and allows quick adaptations to new tasks, to perform survival analysis and pathway enrichment on pan-cancer datasets. In recent machine learning research, meta-learning has been effectively used for knowledge transfer among multiple related datasets. RESULTS: We use meta-learning with Cox hazard loss to show that the integration of TCGA pan-cancer data increases the performance of survival analysis. We also apply advanced model interpretability method called DeepLIFT (Deep Learning Important FeaTures) to show different sets of enriched pathways for multi-omics and transcriptomics data. Our results show that multi-omics cancer survival analysis enhances performance compared with using transcriptomics or clinical data alone. Additionally, we show a correlation between variable importance assignment from DeepLIFT and gene coenrichment, suggesting that genes with higher and similar contribution scores are more likely to be enriched together in the same enrichment sets. AVAILABILITY AND IMPLEMENTATION: https://github.com/berkuva/TCGA-omics-integration.

Perturb-tracing enables high-content screening of multiscale 3D genome regulators.

Three-dimensional (3D) genome organization becomes altered during development, aging, and disease1-23, but the factors regulating chromatin topology are incompletely understood and currently no technology can efficiently screen for new regulators of multiscale chromatin organization. Here, we developed an image-based high-content screening platform (Perturb-tracing) that combines pooled CRISPR screen, a new cellular barcode readout method (BARC-FISH), and chromatin tracing. We performed a loss-of-function screen in human cells, and visualized alterations to their genome organization from 13,000 imaging target-perturbation combinations, alongside perturbation-paired barcode readout in the same single cells. Using 1.4 million 3D positions along chromosome traces, we discovered tens of new regulators of chromatin folding at different length scales, ranging from chromatin domains and compartments to chromosome territory. A subset of the regulators exhibited 3D genome effects associated with loop-extrusion and A-B compartmentalization mechanisms, while others were largely unrelated to these known 3D genome mechanisms. We found that the ATP-dependent helicase CHD7, the loss of which causes the congenital neural crest syndrome CHARGE24 and a chromatin remodeler previously shown to promote local chromatin openness25-27, counter-intuitively compacts chromatin over long range in different genomic contexts and cell backgrounds including neural crest cells, and globally represses gene expression. The DNA compaction effect of CHD7 is independent of its chromatin remodeling activity and does not require other protein partners. Finally, we identified new regulators of nuclear architectures and found a functional link between chromatin compaction and nuclear shape. Altogether, our method enables scalable, high-content identification of chromatin and nuclear topology regulators that will stimulate new insights into the 3D genome functions, such as global gene and nuclear regulation, in health and disease.

Intrinsic bias estimation for improved analysis of bulk and single-cell chromatin accessibility profiles using SELMA.

Genome-wide profiling of chromatin accessibility by DNase-seq or ATAC-seq has been widely used to identify regulatory DNA elements and transcription factor binding sites. However, enzymatic DNA cleavage exhibits intrinsic sequence biases that confound chromatin accessibility profiling data analysis. Existing computational tools are limited in their ability to account for such intrinsic biases and not designed for analyzing single-cell data. Here, we present Simplex Encoded Linear Model for Accessible Chromatin (SELMA), a computational method for systematic estimation of intrinsic cleavage biases from genomic chromatin accessibility profiling data. We demonstrate that SELMA yields accurate and robust bias estimation from both bulk and single-cell DNase-seq and ATAC-seq data. SELMA can utilize internal mitochondrial DNA data to improve bias estimation. We show that transcription factor binding inference from DNase footprints can be improved by incorporating estimated biases using SELMA. Furthermore, we show strong effects of intrinsic biases in single-cell ATAC-seq data, and develop the first single-cell ATAC-seq intrinsic bias correction model to improve cell clustering. SELMA can enhance the performance of existing bioinformatics tools and improve the analysis of both bulk and single-cell chromatin accessibility sequencing data.

Nonparametric single-cell multiomic characterization of trio relationships between transcription factors, target genes, and cis-regulatory regions.

The epigenetic control of gene expression is highly cell-type and context specific. Yet, despite its complexity, gene regulatory logic can be broken down into modular components consisting of a transcription factor (TF) activating or repressing the target gene expression through its binding to a cis-regulatory region. We propose a nonparametric approach, TRIPOD, to detect and characterize the three-way relationships between a TF, its target gene, and the accessibility of the TF's binding site using single-cell RNA and ATAC multiomic data. We apply TRIPOD to interrogate the cell-type-specific regulatory logic in peripheral blood mononuclear cells and contrast our results to detections from enhancer databases, cis-eQTL studies, ChIP-seq experiments, and TF knockdown/knockout studies. We then apply TRIPOD to mouse embryonic brain data and identify regulatory relationships, validated by ChIP-seq and PLAC-seq. Finally, we demonstrate TRIPOD on the SHARE-seq data of differentiating mouse hair follicle cells and identify lineage-specific regulation supported by histone marks and super-enhancer annotations. A record of this paper's transparent peer review process is included in the supplemental information.

Extracting physical characteristics of higher-order chromatin structures from 3D image data.

Higher-order chromatin structures have functional impacts on gene regulation and cell identity determination. Using high-throughput sequencing (HTS)-based methods like Hi-C, active or inactive compartments and open or closed topologically associating domain (TAD) structures can be identified on a cell population level. Recently developed high-resolution three-dimensional (3D) molecular imaging techniques such as 3D electron microscopy with in situ hybridization (3D-EMSIH) and 3D structured illumination microscopy (3D-SIM) enable direct detection of physical representations of chromatin structures in a single cell. However, computational analysis of 3D image data with explainability and interpretability on functional characteristics of chromatin structures is still challenging. We developed Extracting Physical-Characteristics from Images of Chromatin Structures (EPICS), a machine-learning based computational method for processing high-resolution chromatin 3D image data. Using EPICS on images produced by 3D-EMISH or 3D-SIM techniques, we generated more direct 3D representations of higher-order chromatin structures, identified major chromatin domains, and determined the open or closed status of each domain. We identified several high-contributing features from the model as the major physical characteristics that define the open or closed chromatin domains, demonstrating the explainability and interpretability of EPICS. EPICS can be applied to the analysis of other high-resolution 3D molecular imaging data for spatial genomics studies. The R and Python codes of EPICS are available at https://github.com/zang-lab/epics.

Single-cell chromatin profiling of the primitive gut tube reveals regulatory dynamics underlying lineage fate decisions.

Development of the gastrointestinal system occurs after gut tube closure, guided by spatial and temporal control of gene expression. However, it remains unclear what forces regulate these spatiotemporal gene expression patterns. Here we perform single-cell chromatin profiling of the primitive gut tube to reveal organ-specific chromatin patterns that reflect the anatomical patterns of distinct organs. We generate a comprehensive map of epigenomic changes throughout gut development, demonstrating that dynamic chromatin accessibility patterns associate with lineage-specific transcription factor binding events to regulate organ-specific gene expression. Additionally, we show that loss of Sox2 and Cdx2, foregut and hindgut lineage-specific transcription factors, respectively, leads to fate shifts in epigenomic patterns, linking transcription factor binding, chromatin accessibility, and lineage fate decisions in gut development. Notably, abnormal expression of Sox2 in the pancreas and intestine impairs lineage fate decisions in both development and adult homeostasis. Together, our findings define the chromatin and transcriptional mechanisms of organ identity and lineage plasticity in development and adult homeostasis.

Single-nucleus chromatin accessibility profiling highlights regulatory mechanisms of coronary artery disease risk.

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across cell types. Genome-wide association studies have identified over 200 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis-regulatory elements. Here, we applied single-nucleus assay for transposase-accessible chromatin with sequencing to profile 28,316 nuclei across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped ~320,000 accessible sites across all cells, identified cell-type-specific elements and transcription factors, and prioritized functional CAD risk variants. We identified elements in smooth muscle cell transition states (for example, fibromyocytes) and functional variants predicted to alter smooth muscle cell- and macrophage-specific regulation of MRAS (3q22) and LIPA (10q23), respectively. We further nominated key driver transcription factors such as PRDM16 and TBX2. Together, this single-nucleus atlas provides a critical step towards interpreting regulatory mechanisms across the continuum of CAD risk.